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A  Multiphasic  Investigation  of  an  Argon  Plasma 
in  a  rf  Discharge  and  Sonic  Afterglow 

Louis  N.  Medgyesi-Mitschang  and  Ray  A.  Hefferlin* 
McDonnell  Douglas  Research  Laboratories 
McDonnell  Douglas  Corporation,  St.  Louis,  Missouri  63166 


Abstract 


Microwave  interferometric  and  radiometric,  optica!  spectrometric  and  thermo¬ 
couple  measurements  of  the  electron  density  and  the  electron,  the  excitation,  and 
the  gas  temperatures  have  been  made  in  the  active  discharge  and  sonic  afterglow  of 
an  argon  plasma  generated  by  a  rf  generator..  The  optical  measurements  were 
made  in  the  4000-5000  A  region  for  discharge  pressures  of  60-650  Torr.  Excita¬ 
tion  temperatures  in  the  discharge  ranged  between  ~7000-l0,000  K.  The  electron 
temperatures  from  the  radiometer  measurements  were  compared  with  the  excita¬ 
tion  temperatures.  Thermocouple  measurements,  gas  dynamic  calculations  and 
rotational  spectra  from  traces  of  N2  injected  in  the  argon  indicated  that  the  gas 
temperatures  in  the  discharge  and  the  afterglow  were  29004400  K  and  1900-3300 
K,  respectively,  indicating  a  two-temperature  state  in  both  regions.  Results  of  this 
investigation  indicated  local  thermodynamic  equilibrium  (LTE)  for  all  excited  levels. 
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1 1ntroduction  and  background 


The  processes  which  govern  the  production  and 
loss  of  electrons  together  with  reactions  between 
ions  and  neutrals  play  an  important  role  in  deter¬ 
mining  the  plasma  environment  of  a  vehicle  during 
planetary  reentry.  In  many  instances  it  would  be 
desirable  to  be  role  to  control  the  composition  of 
this  plasma.  To  accomplish  this  objective,  it  is 
necessary  to  know  what  are  the  dominant  reaction 
processes  and  reaction  rates  at  ter  peratures  en¬ 
countered  during  reentry  (20004000  K).  Also  of 
interest  are  the  chemical  reactions  which  occur 
when  target  gas  molecules  interact  with  a  plasma 
in  this  high  temperature  range. 

Reactions  associated  with  plasma,  at  elevated 
temperatures  (>  1000  K)  exhibit  significantly 
different  behavior  than  the  analogous  reactions 
occurring  at  room  temperature.  The  specific  reac¬ 
tion  pathways,  and  rate  constants  of  these  reac¬ 
tions  are  of  special  interest  at  these  temperatures/'2 

An  experimental  system  has  been  developed  at 
McDonnell  Douglas  Research  Laboratories  tnat 
permits  investigation  of  high  temperature  plasma 
reactions  over  a  broad  range  of  plasma  conditions.2 
At  present,  only  electron  attachment  reactions  are 
being  investigated  because  of  the  wide  interest  in 
electrophilic  molecules.  A  sonic  afterglow  serves 
as  a  chemical  reaction  region  for  interactions  be¬ 
tween  the  neutral  electrophilic  molecules  and  the 
plasma.  The  state  of  the  plasma  is  monitored  by  a 
variety  of  diagnostic  techniques  in  both  the  dis¬ 
charge  and  afterglow  regions.  The  results  of  opti¬ 
cal  spectrometer,  microwave  radiometer,  micro¬ 


wave  interferometer  measurements,  as  well  as 
supplementary  tncrmocouple  data  and  gas  temper¬ 
ature  calculation,  presented  in  this  paper  allow 
accurate  characterization  of  the  conditions  in  the 
discharge  and  afterglow  regions.  Knowledge  of 
these  conditions  is  essential  for  the  interpretation 
of  the  chemical  kinetics  involved  in  the  reaction 
between  the  plasma  and  the  injected  electrophilic 
molecules. 

In  the  past,  many  techniques  have  been  developed 
to  determine  the  important  plasma  parameters  such 
as  the  electron  temperature,  excitation  temperature, 
and  electron  density.  Most  prominent  of  these 
have  been  Langmuir  probes  and  optical  spectro- 
metric  methods.  Spectrometric,  radiometric,  and 
interferometric  techniques  haw  the  advantage  of 
allowing  plasma  diagnosis  without  introducing 
perturbations  into  the  gas.  Early  experiments  by 
Knot* ,  and  Easley  and  Munford-5  have  demon¬ 
strated  that  the  radiometrically  determined  electron 
temperatures  compared  favorably  with  those  ob¬ 
tained  from  Langmuir  probe  measurements  in  low 
pressure  Hg-Ar  discharges.  More  recently,  McGregor 
and  Brewer6  have  found  that  in  certain  cases  the 
spectrometrically  measured  excitation  temperatures 
correlate  with  the  electron  temperatures  measured 
with  Langmuir  probes  in  an  argon  plasma  jet. 

Similar  correlation  of  the  radiometric  temperature 
data  with  spectrometric  results  was  found  by  Dellis7 
and  Harding  et  al.* 

In  this  report  the  excitation  temperature  meas¬ 
ured  by  a  spectrometer  and  the  electron  tempera¬ 
ture  obtained  from  a  radiometer  are  compared. 
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INTRODUCTION  AND  BACKGROUND 


Electron  densities  obtained  from  spectrometric 
measurements  in  the  sonic  afte  glow  channel  are 
compared  with  interferometrically  measured  values 
extrapolated  into  the  channel  from  the  free  jet. 

The  calculated  gas  (heavy  particle)  temperatures 
obtained  from  gas  dynamic  considerations  both  in 
the  channel  and  the  discharge  chamber  and  from 
molecular  rotational  lines  in  the  channel,  are  com¬ 
pared  with  thermocouple  data  obtained  in  the  after¬ 
glow  channel  for  various  levels  of  input  rf  power. 
These  measurements  are  used  to  examine  the 
existence  of  local  thermodynamic  equilibrium 
(LTE)  in  the  active  discharge  region  and  the  sonic 
afterglow. 

LTE  is  expected  to  hold  if  collisional  processes 
involving  a  Maxwellian  electron  distribution  dom¬ 
inate  the  rate  equations  governing  the  population 
of  a  given  quantum  level.  A  “critical  level”  may 
be  defined  at  which  the  collisional  excitation  rates 
are  comparable  to  or  exceed  the  radiative  decay 


rates.  The  (electron)  temperature  associated  with 
the  energy  distribution  of  free  electrons  can  be 
assumed  to  be  in  equilibriun  with  the  (excitation) 
temperature  of  states  whose  levels  are  above  the 
critical  le”'*'. 

If  LTE  is  established  above  a  certain  level, 
the  population  in  the  various  levels  above  tnis  level 
can  be  related  with  each  other  via  the  Boltzmann 
distribution..  In  addition,  if  it  can  be  shown  that 
the  phsma  is  optically  thick,  reabsorption  of  radia¬ 
tion  i  an  occur,  resulting  in  an  effective  reduction  of 
the  radiative  decay  rates  so  that  LTE  can  exist  at 
some  lower  level  than  tnat  at  which  radiation  and 
collisions!  processes  are  co  mparable.  However,  as 
will  be  noted  later,  the  plasma  in  both  the  active 
discharge  and  the  sonic  afteiplow  (reaction  channel) 
regions  is  optically  thin.  Conditions  mitigating 
against  unqualified  LTE  are  time  variations  in  the 
discharge,  spat.al  gradients  in  the  electron  density, 
and  strong  electromagnetic  fields. 
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2  Experimental  apparatus 


The  apparatus  associated  with  this  experiment  is  scribed  here.  The  experimental  apparatus  germane 

shown  in  Fig.  1 .  It  has  been  described  in  consider-  to  this  discussion  consists  of  the  plasma  source 
able  detail  elsewhere/  Only  those  elements  of  the  followed  by  tlr  reaction  channel,  the  spectrometer 

system  which  have  a  direct  bearing  upon  the  spec-  and  ancillary  optical  system,  and  the  microwave 

f  r'>metric  and  microwave  measurements  will  be  de-  radiomtier  and  microwave  interferometer. 


roughing  diffusion 

pump  pump 


Fig.  1  Experimental  apparatus  for  high  temperature  plasma  reaction  studies 
(Mullen  et  a!..  Rev.  Sci,  Instr.  41, 1746  (1970)) 
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EXPERIMENTAL  APPARATUS 


For  the  source,  an  electrodeless  induction  argon 
plasma  similar  to  that  of  Reed^and  Eckert  el  al.70, 
is  used.  This  source  is  free  of  electrode  contamina¬ 
tion  making  it  particularly  suitable  for  high  temper¬ 
ature  chemistry  studies.  The  plasma  generator  is  a 
100  kW,  class  C,  rf  oscillator  operating  at  540  kHz. 
The  power  input  to  the  rf  oscillator  can  be  varied 
continuously  from  0  to  100  kW.  For  the  measure¬ 
ments  described  in  this  paper  input  rf  power  levels 
of  42, 48,  54,  and  60  kW  were  chosen.  The  following 
variables  are  monitored:  oscillator  plate,  grid,  tank 
voltages  and  currents,  input  power  per  phase,  grid 
coupling,  feedback  control  and  rf  frequency.  The 
temperature  rise  in  the  rf  coil  cooling  water,  the 
chamber  pressure  anJ  the  input  mass  flow  rate  are 
also  monitored.  Any  generator  operating  condition 
can  be  reset  to  ±  5%  duplicating  previously  obtained 
conditions. 


The  discharge  tube  is  made  of  1  mm  thick  quartz, 
80  mm  in  diam  160  mm  long,  and  is  sandwiched  be¬ 
tween  0-rin",ed,  water-cooled  copper  end  plates. 

The  relative  physical  dimensions  of  the  discharge 
chamber  itself  are  indicated  in  Fig.  2.  The  dis¬ 
charge  gas  inlets  pass  through  one  of  the  copper 
end  plates.  Typical  argon  mass  flow  rates  are  0.10 
to  0.36  g/sec  as  shown  in  Fig.  3.  The  rf  power  is 
coupled  through  two  eight-turn  counterwound 
coils  of  3  mm  o.d.  copper  tubing  that  surround  the 
plasma  tube.  The  coils  are  center  fed  for  minimum 
rf  plasma  potential  at  each  end  of  the  discharge.7 7 
A  1 .5  mm  clearance  between  the  discharge  tube  ana 
a  Teflon  water  jacket  allows  adequate  coolant  water 
flow.  The  partially  ionized  argon  is  passed  from 
the  discharge  chamber  into  a  4  mm  diam  sonic  re¬ 
action  channel,  and  then  expanded  as  a  free  jet  into 
a  large  vacuum  tank  at  0.2  Torr  pressure. 


To  reaction 
channel 
(not  shown) 


"Oppositely  wound,  8  turns 


0  5  10  15 

I _ _ _ I _ _ _ 1 - - — J 

Scale  ■  cm 


Fig.  2  Discharge  chamber  and  arrangement  of  rf  excitation  coils 
(Mullen  et  al..  Rev.  Sci.  Instr.  41, 1746  (1970)) 
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0 1 _ I _ I - 1 - 1 

0  0.1  0.2  0.3  0.4 

Mass  flow  rate,  M  (g/sec) 


Fig.  3  Measured  mass  flow  rates  in  discharge 
chamber 


2.1  Optical  instrumentation 


Spectroscopic  determination  of  the  argon  excita¬ 
tion  temperature  and  the  electron  density  in  the 
discharge  chamber  and  the  reaction  channel  was 
made  using  the  optical  arrangement  outlined  in 
Fig.  4.  A  1  m  (Jarrell-Ash  Model  78-466) 
Czerny-Turner  scanning  spectrometer  was  used  in 
these  experiments.  The  plasma  radiation  from  the 
discharge  chamber  was  transferred  to  the  spectro¬ 
meter  using  a  0.32  cm  diam  by  14.1  cm  long  fiber 
bundle,  with  plane-polished  ends.  The  bundle  was 
mounted  in  the  radiometer  port  with  its  input  end 


flush  against  the  quartz  window.  Then,  a  lens  was 
used  to  form  an  image  of  the  end  of  the  bundle  on 
the  spectrometer  slit.  A  similar  arrangement  with 
an  aluminized  quartz  rod  was  tried  and  found  to  be 
less  satisfactory  due  to  the  fragility  of  the  aluminum 
coating.  A  photomultiplier  tube  was  used  w:*h  the 
resulting  output  displayed  by  a  strip  chart  recorder 
(modified  Honeywell  Electronic  1 6). 

A  plane  window  standard  35A  tungsten  lamp 
was  used  for  calibration  of  the  optical  system.  The 
entire  measurement  system  from  the  fiber  bundle 
onward,  was  removed  from  the  discharge  chamber 
and  the  light  from  the  lamp  was  transferred  to  the 
input  of  the  fiber  bundle.  Spectrum  scans  were  ob¬ 
tained  in  the  4000-5000  A  region.  A  filter  was  used 
to  eliminate  second-order  features. 

Spectroscopic  measurements  in  the  discharge 
were  made  for  a  range  of  discharge  pressures  of 
60-650  Torr  and  input  power  levels  of  42-60  kW. 

The  spectrum  scans  were  obtained  over  the  entire 
pressure  range  by  adjusting  the  mass  flow  rate  to 
the  discharge  chamber,  holding  the  input  power 
constant.  Usually  the  spectrum  scans  were  made 
for  increasing  chamber  pressures,  followed  by  scans 
retracing  the  pressure  range..  In  this  way  the  repeat¬ 
ability  of  the  data  could  be  checked.  During  data 
acquisition,  the  rf  generator  was  operated  uninter¬ 
rupted  for  long  periods  (often  exceeding  one  hour). 
During  these  long  periods  of  operation  the  genera¬ 
tor  was  found  to  be  stable  and  showed  no  percep¬ 
tible  drift  in  its  operating  characteristics.  The 
special  anti-sputtering  design  of  the  rf  coil  and  the 
optimized  configuration  of  the  discharge  chamber 
prevented  contamination  of  the  plasma. 

For  the  optical  measurements  in  the  reaction 
channel,  the  fiber  bundle  light  pipes  were  integrally 
bonded  to  the  reaction  channel..  Their  mounting 
was  flush  with  the  inside  wali  of  the  channel,  causing 
no  perturbation  in  the  plasma  stream.  Light  inten¬ 
sities  objerved  were  substantially  weaker  than  those 
from  the  discharge  chamber  and  care  was  taken  to 
align  the  spectrometer  via  a  small  laser.  Again, 
spectrum  scans  were  made  for  the  previously  noted 
range  of  rf  power  and  discharge  chamber  pressures. 
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EXPERIMENTAL  APPARATUS 


3  high) 


2.2  Radiometer  and  associated 
instrumentation 

For  the  radiometric  measurements  a  modified  Dicke 
X-band  microwave  radiqmeter  operating  at  8.5  GHz 
center  frequency  was  used,12'15  The  basic  radiom¬ 
eter  circuit  is  shown  in  block  diagram  form  in  Fig.  5. 
In  operation  the  radiometer  was  electronically 
switched  between  the  plasma  and  a  standard  wave¬ 
guide  noise  source  consisting  of  an  argon-filled  tube 
mounted  obliquely  to  the  E-plane  in  a  rectangular 
waveguide.  The  standard  source  presents  a  low 
VSWR,  less  than  1,2,  with  or  without  an  argon 


plasma  present  in  the  above  tube.  Spurious  electro¬ 
magnetic  modes  due  to  the  presence  of  the  gas  tube 
are  suppressed  within  the  waveguide.  With  the  gas 
tube  on,  the  white  noise  spectrum  is  generated 
over  the  bandwidth  of  the  waveguide.  For  radiom¬ 
etric  measurements,  the  gas  tube  was  operated  in 
the  cw  mode.  The  standard  noise  source,  together 
with  a  50  dB  precision  attenuator  formed  a  cali¬ 
brated  blackbody  temperature  standard  with  equiv¬ 
alent  temperature  >  10,000  K.  Due  to  the  low 
plasma  emissivity,  the  effective  thermal  radiation 
from  the  plasma  corresponded  to  a  substantially 
lower  temperature. 


6 


MCDONNELL  DOUGLAS 


Report  MDC  <2046 1  '■  September  19?  1 


EXPERIMENTAL  APPARATUS 


The  signal  from  the  plasma  and  the  standard  are 
alternately  gated  into  the  radiometer  by  a  crystal 
switch.  The  two  signals  ar?  detected,  amplified, 
and  processed  by  passing  them  through  a  balanced 
low  noise  mixer  i.f.-pre-amp  combination  followed 
by  a  high  gain  i.f.  amplifier  and  a  coherent  detector 
synchronized  with  the  crystal  switch,  as  shown  in 
Fig.  5.  The  output  is  detected  on  a  precision  micro¬ 
voltmeter. 

To  assure  symmetry  in  operation,  test  signals 
were  fed  to  ports  A  and  B  of  the  crystal  switch, 
alternately  switched  and  compared.  Slight  discrep¬ 
ancies  between  the  two  ports  were  corrected  by  ad¬ 
justing  the  E-H  tuning  plungers  in  the  switch.  Bal¬ 
anced  operation  was  verified  by  connecting  identi¬ 
cal  noise  sources  to  the  input  ports  of  the  switch 
and  observing  a  null  on  the  radiometer  output. 

The  radiometer  receiving  antenna  is  a  contoured 
open-ended  waveguide,  japped  t>v  a  similarly  con¬ 
toured  1.5  mm  thick  quartz  window.  The  wave¬ 
guide  is  oriented  so  that  the  E-plane  is  parallel  to 
the  axis  of  the  discharge  chamber,  resulting  in  the 
minimum  mismatch.  The  alternate  polarization 
proved  less  satisfactory  for  coupling.  Contouring 
the  waveguide  to  the  curvature  of  the  discharge 
tube  minimized  the  necessary  water  gap  between 
the  antenna  and  the  discharge  jacket,  thus  limiting 
the  absorption  in  the  water  film  to  1  dB.  The  above 
optimum  arrangement  satisfied  the  dual  requirements 
of  close  coupling  of  the  antenna  to  the  plasma  and 
proper  cooling  of  both  antenna  window  and  dis¬ 
charge  tube. 

Initially  with  the  discharge  off,  the  radiometer 
antenna  with  an  E-H  tuner  was  connected  to  the 
discharge  chamber  using  a  coherent  test  signal  at 
the  center  frequency  of  the  radiometer.  An  isolator 
was  inserted  to  maintain  a  low  VSWR  in  the  antenna 
leg  with  changing  plasma  conditions  and  to  prevent 
any  local  oscillator  leakage.  Upon  switching  to  the 
radiometer  mode,  a  signal  corre^oonding  to  a  300  K 
blackbody  source  was  measured.  With  the  plasma 
on,  the  thermal  radiation  from  the  plasma  was 
measured  radiometrically  at  a  generator  input  power 
of  54  kW,  representing  the  mid-range  of  input  power.. 


In  these  measurements  the  pressure  in  the  discharge 
chamber  was  varied  by  changing  the  mass  flow  rate 
in  the  chamber.  The  rf  power  was  maintained  con¬ 
stant.  The  readings  were  obtained  by  setting  the 
precision  attenuator  for  a  null  output.  The  average 
of  three  consecutive  readings  was  taken  at  each 
pressure  setting  and  the  result  recorded.  After 
sweeping  through  the  pressure  range,  as  with  the 
optical  measurements,  the  pressure  ran,  ,e  was 
retraced  and  the  data  were  compared  with  previ¬ 
ously  obtained  values  for  repeatability.  Suitable 
cooling  of  the  discharge  chamber  allowed  all  the 
microwave  components  in  close  proximity  with  the 
discharge  to  remain  at  a  constant  ambient  tempera¬ 
ture  during  the  experiment.  After  a  sweep  plus  re¬ 
trace  of  the  pressure  range,  the  discharge  was  extin¬ 
guished  and  the  radiometer  nulled  with  the  preci¬ 
sion  attenuator..  The  attenuator  setting  was  com¬ 
pared  with  the  initial  value  at  the  beginning  of  the 
experiment,  thereby  serving  as  an  index  of  drift 
stability.  Because  of  the  good  drift  compensation 
charactesistics  of  the  mixer  and  i.f.,  stages,  negligi¬ 
ble  drut  was  noted  when  the  microwave  circuit 
was  properly  tuned. 

To  obtain  an  electron  temperature  from  the  radi¬ 
ometer  measurements,  the  plasma  emissivity  was 
determined.:  A  plane  wave  analysis  compatible  with 
the  physical  dimensions  of  the  radiometer  antenna 
and  the  diameter  of  the  plasma  discharge  (8  cm  in 
diam)  was  used.  Incoherence  in  the  power  reflec¬ 
tance  inside  the  plasma  was  assumed.  Separate 
reflectance  and  transmission  measurements  were 
made  in  the  discharge  at  various  chamber  pressures. 
The  plasma  induced  attenuation  was  shown  to  ex¬ 
ceed  50  dB.  Using  these  considerations,  the  expres¬ 
sion  for  the  plasma  emissivity  was  computed  as  a 
function  of  the  total  pressure  in  the  discharge,  the 
gas  temperature,  the  measured  mass  flow  in  the 
chamber,  a^  1  flie  electron  density..  The  analysis 
included  ,:ne  effect  of  both  Coulomb  and  electron- 
neutral  interactions.  The  main  points  of  this 
analysis  are  summarized  in  Appendix  \. 
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Fig.  5  Radiometer  schematic 
(Mullen  et  al..  Rev.  Sci.  Ins tr.  41, 1746  (1970)) 


EXPERIMENTAL  APPARATUS 


2.3  Microwave  interferometer 

Electron  densities  were  measured  in  the  free  jet  - 
region,  3.S8  cm  downstream  from  the  exit  plane  of 
the  reaction  channel  using  a  microwave  interferom¬ 
eter  as  shown  in  Fig.  6.  Open  ended  wave  guides 
(E  plane  polarized  perpendicular  to  the  channel 
axis)  used  for  antennas  are  located  so  that  the 
plasma  is  just  outside  of  the  antenna  near  field.  To 
prevent  reflections  in  the  plasma  leg  of  the  interfer¬ 
ometer  from  entering  the  circuit,  an  isolator  is  used 
behind  the  transmitting  antenna.  Similarly,  another 
isolator  is  placed  behind  the  receiving  antenna  so 
that  the  crystal  detector  can  be  tuned  to  the  inci¬ 
dent  signal  only,  avoiding  reactance  tuning  due  to 
the  plasma  path  length.  Precision  phase  shifters  and 
attenuators  are  used  throughout  the  circuit  to  ob¬ 
tain  maximum  accuracy. 

With  the  plasma  off,  the  interferometer  is  null 
adjusted  using  the  precision  attenuators  and  phase 


shifters.  With  the  plasma  on,  the  path  propagation 
characteristics  between  the  antennas  is  altered  and 
the  interferometer  is  renulled.  The  plasma  by-pass 
arm  in  the  circuit  permitted  tht  interferometer  to 
be  nulled  during  the  experiment,  thereby  checking 
drift.  Because  the  plasma  is  underdose  at  the  exit 
plane  where  the  interferometer  is  sitiiated  (Fig.  1), 
the  electron-ion  collision  frequency  v  is  small  com¬ 
pared  to  the  interferometer  signai  frequency  f  (Le., 
<!  10~3).  No  appreciable  plasma  attenuation  is 
observed  and  the  electron  density  may  be  deter¬ 
mined  from  phase  shift  measurements  alone,  using 
a  collisionless  plasma  model  with  a  parabolic  profile. 

Assuming  that  isentropic  expansion  occurs  past 
the  exit  plane  of  the  reaction  channel,  the  electron 
density  was  shown  to  be  approximately  a  factor  of 
354  higiier  in  the  channel  than  t^at  measured  3.58 
cm  further  downstream.-*  The  electron  densities 
were  measured  over  the  range  of  conditions  paral¬ 
leling  those  for  the  optical  measi  rements. 


Plasma 

Bypass 


Fig.  6  Ka-band  microwave  interferometer  block  diagram 
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3  Summary  of  measurements 


3.1  Optical  measurements 


For  the  chamber  pressures  investigated  the  optical 
data  contained  no  visible  Aril  lines  in  the  spectra. 
The  Arl  were  in  l.n  case  wider  than  the  instrument 
resolution  width  ( 1 .2  A).  This  observation  allowed 
placing  an  order  of  magnitude  upper  limit  of  1 0  * 6 
e/cm  on  the  '•'sctron  density.  For  excitation  tem¬ 
perature  determination,  the  following  Arl  lines 
were  selected:  4158.6  A,  4164.2  A,  4181.9A, 

4259.4  A,  4266.3  A,  4272.2  A,  4300. 1  A,  4345.2  A, 
and  4510.7  A.  Boltzmann  plots  were  obtained 
f.om  the  measured  line  intensities.  From  these 
the  excitation  temperatures  were  obtained  by 
least  squares  fitting  to  the  Boltzmann  relationship 
describing  the  intensity  of  a  spectral  line, 


.  1X 

,ogj0  f  -'o?  10 

8  nm 


5040 

T 

‘exc 


where  h  = 
g  = 


o 

X  = 
I  = 


^n  “ 
^exc  = 


Planck’s  constant, 
statistical  weight, 
emission  oscillator  strength, 
atom  partition  function, 
total  number  of  atoms  per  cm3, 
wavelength  of  emitted  ratiation, 
intensity  of  atomic  spectral  line 
(ergs  sec-'  ster-'  cm  -2), 
excitation  energy  of  level  n  (eV), 
and 

excitation  temperature. 


In  usini  this  relationship  two  assemptions  were 
implicit:  1 )  The  plasma  war  assumed  to  be  optically 
thin  in  the  wavelength  region  scanned  so  that  self¬ 
absorption  was  absent  and  2)  ti:°  excited  states 
were  assumed  to  possess  a  Boltzmann  distribution. 

The  plasma  was  shown  to  be  optically  thin  and 
thus  no  self-absorption  was  present  since  the  meas¬ 
ured  spectral  radiances  of  the  Arl  lines  were  not 
comparable  to  the  values  predicted  by  the  Planck 
function.  For  example,  the  peak  measured  spectral 
radiances  were  It)'  '  -lo'3  ergs  sec-'  ster-'  cm-3. 
The  Planck  function  at  the  same  wavelengths  pre¬ 
dicts  ~  2  x  10' 3  ergs  sec-'  ster-'  cm-3  for 
T  =  9000  K  in  the  discharge  and  ~1  x  lo'^  ergs 
sec-'  ster-'  cm-3  for  T  =  5000  K  in  the  reaction 
channel. 

The  excitation  temperatures  in  the  discharge 
chamber  and  the  reaction  channel  are  summarized 
in  Tables  1  and  2.  Figure  7  shows  a  typical  Boltzmann 
plot  for  the  data  obtained  in  the  reaction  channel. 

The  electron  density  was  calculated  from  the 
continuum  radiation  via  the  recently  corrected 
Kramers-Unsold  equation/6  This  formulation 
states  that  the  spectral  radiance  of  the  continuum 
of  an  LTE  argon  plasma  at  an  electron  temperature 
Te  (and  with  a  calculable  electron  density,  Ne)  is 

ev  =  (5.443)  x  10~39  Ne2  T“,/2 
ergs  sec-'  ster-'  cm-3  Hz. 
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Tabii*  1  Representative  spectrometric  results  for  discharge  chamber 


Chamber 

Power 

Taxc 

N. 

Chamber 

Power 

1 

(Torr) 

(kW) 

IK) 

(•/cm3) 

(Torr) 

(kW) 

(K) 

(e/cm3) 

75 

42 

9924 

8.5  xIO15 

63 

54 

10218 

1.17  xIO16 

100 

42 

9777 

8.68  xIO15 

100 

54 

9351 

1.17  x  1016 

8.40  x  1015 

127 

54 

9613 

1.16  xIO16 

125 

42 

8778 

150 

54 

9435 

1.02  x  1016 

8.22  x  1015 

180 

54 

9154 

1.17  x  1016 

150 

42 

8062 

•: 

54 

54 

7684 

7592 

0.925  xIO16 
1.14  x  1016 

7.45  x  1015 

1.18  x  1016 

175 

42 

6446 

305 

54 

200 

42 

8234 

7.7  xIO15 

353 

54 

8172 

1.34  x  1016 

250 

42 

8555 

8.63  x  1015 
9.74  x  1015 

417 

54 

9253 

1.43  xIO16 
1.46  x  1016 

■ 

42 

7398 

475 

54 

8223 

42 

8729 

1.0  xIO16 

535 

54 

7352 

1.48  x  1016 

425 

42 

8631 

1.23  x  1016 

644 

54 

8722 

1.51  x  1016 

1.1  x  101 

100 

60 

9729 

1.19  x  1016 
1.07  x  1016 

75 

48 

13709 

150 

8823 

100 

48 

9853 

1.04  x  1016 
8.9  x  1015 

200 

60 

6719 

1.1  x  1016 

1.22  x  1016 

150 

48 

6898 

250 

60 

6064 

1.4  x  1016 

1.27  xIO16 

48 

8389 

300 

60 

7914 

250 

48 

7884 

1.0  xIO16 

300 

48 

7286 

1.09  x  1016 

350 

48 

6723 

1.18  x  1016 

425 

48 

8612 

1.33  x  1016 

525 

48 

7383 

1.37  x  1016 
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Table  2  Representative  spectrometric  results  for  reaction  channel 


E  (eV) 

Y 

14.670 

-9.008 

14.620 

• -8.055 

-8.00 

14  620 

-8.942 

14.510 

-8.884 

14.470 

-8.726 

_0> 

c 

14  470 

-8.796 

o 

14  460 

-8.773 

c 

i 

-8.50 

14.460 

-8  796 

14  440 

-8.789 

2 

>• 

O 

= 4768  K 

CO 


-9.00 


o 

i 


-9.50 


Power 

Che  mber  pressure 


-  argon 
=  60  kW 
=  1 59  T  orr 


14.4 


14.5  14.6 

S  <eV) 


14.7 


Fig.  7  Boltzinann  plot  for  excitation  temperatures 
in  reaction  channel 


The  empirical  constant  in  this  equation  is  taken 
from  Rivhter/  7  and  Olsen7*.  The  continuum  meas¬ 
urements,  taken  at  X  =  4500  A,  were  uninverted, 
that  is,  the  source  was  not  traversed.  Data  taken  for 
different  slit  widths  agreed.  The  oscillator  strengths 
used  in  the  data  reduction  are  given  in  the  recent 
NBS  critical  compilation  of  measured  oscillator 
strengths.7  9 

The  electron  densities  in  the  discharge  chamber, 
using  the  Kramers- Unsold  equation,  are  given  in 
Fig.  8.  These  values  assume  the  effective  depth  of 
the  plasma  to  be  2  cm.  The  measured  densities 
vary  by  a  factor  of  two  (0.70  to  1.4  x  10* 6  e/cm^ ) 
over  the  pressure  range  from  40  to  540  Torr  at  the 
lower  rf  power  levels  (42  and  48  kW),  but  remain  a 
relatively  constant  function  of  discharge  pressure 
at  60  kW.  Similar  values  were  measured  by  Hughes 
and  Wooding’0  and  Goldfarb  and  Dresvin’7  for 
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Table  3  Comparison  of  argon  rf  plasma  sources 


rf  power  coupled  into  plasma 
Synopsis  of  data  given  in  previous  figures 
Center  fed  counter  wound  coil 
Some  channel  diameter 


Electron  density  in  discharge  chamber  (e/cmJ) 


SUMMARY  OF  MEASUREMENTS 


somewhat  higher  frequency  discharges  at  atmos¬ 
pheric  pressure.  But  the  present  values  differ  from 
those  by  Eckert,  et  al.10  For  comparison  of  the 
experimental  parameters  and  resultant  oata  of 
these  investigators  see  Table  3. 


Fig.  8  Electron  density  in  discharge  chamber  using 
spectroscopic  data  and  the  Kramers-Unsold 
formulation 


3.2  Other  measurements 

The  interferometer  measurements  show  the  electron 
density  in  the  channel  to  be  a  monotonically  in¬ 
creasing  function  of  the  discharge  pressure  over  the 
range  considered  (Fig.  9).  Similar  behavior  is  noted 
for  the  electron  density  as  a  function  of  the  average 
reaction  channel  temperatures  determined  via 
thermocouple  measurements  for  42,  48,  and  54  kW, 
as  shown  in  Fig.  10. 

Using  the  Kramers-Unsold  equation,  spectrom- 
etric  data  in  the  channel  yield  electron  densities 
that  are  in  good  agreement  with  interferometer 
measurements  extrapolated  upstream  from  the 
point  of  measurement  in  the  free  jet  (Fig.  1  la- 1  Id). 


The  fact  that  the  optical  and  microwave  methods 
yield  values  within  25%  over  a  wide  range  of  plasma 
conditions  is  significant,  considering  that  the  meas¬ 
urement  techniques  are  based  on  completely 
different  principles. 


Fig.  9  Reaction  channel  electron  density  using 
interferometer  measi-ements 

Values  for  the  plasma  emissivity  were  obtained 
by  taking  the  ratio  of  the  effective  blackbody 
radiation  and  excitation  temperatures  at  various 
discharge  pressures.  This  procedure  assumed  that 
the  electron  and  excitation  temperatures  were  com¬ 
parable  due  to  the  efficient  energy  exchange  be¬ 
tween  free  electrons  and  the  upper  bound  electronic 
states.  In  Fig.  1 2,  the  values  for  the  emissivity  are 
compared  with  calculated  ones  obtained  from  a  slab 
plasma  model  (Appendix  A).  In  this  calculation  the 
spectrometrically  determined  values  for  the  excita¬ 
tion  temperatures  were  used  for  the  electron  tem¬ 
peratures.  Agreement  between  the  calculate]  and 
measured  emissivities  (Fig.  1 2)  indicates  that  the 
excitation  and  electron  temperatures  approximate 
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each  other  in  the  discharge.  The  electron  densities 
from  the  emisshrity  calculations  shown  in  Fig.  1 2 
are  compared  with  the  optical  measurements.  The 
computed  values  of  the  electron  density  are  higher 


than  the  measured  values.  Using  the  measured 
electron  density  directly  in  the  emissivity  program 
did  not  significantly  alter  the  overall  emissivity 
results  for  the  range  of  discharge  pressures  considered. 


Fig.  10  Reaction  channel  electron  density  as  a  function  of  channel  gas  temperature.  Electron  density 
measured  with  Ka-band  interferometer.  Gas  temperature  measured  with  a  tungsten-rhenium 
thermocouple  placed  in  the  centerline  of  the  channel 
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Electron  density  in  reaction  channel,  x  1010 


SUMMARY  OF  MEASUREMENTS 


Discharge  chamber  pressure  (Torr) 


FIG.  11  Comparison  of  spectrometer  and  interferometer  data  in  reaction  channel 


Discharge  chamber  pressure  (Torr)  uj 


Fig.  12  Plasma  emissivity  and  electron  density  as  a 
function  of  discharge  pressure  rf  power  -  54  kW 
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4  Discussion 


The  state  of  the  plasma  in  the  discharge  and  the 
sonic  afterglow  (reaction  channel)  regions,  as 
determined  from  the  preceding  measurements 
are  summarized  below.  Additional  discussion  of 
thermocouple  results  as  well  as  supplementary 
optical  and  gas  dynamic  considerations  are  included. 

(a)  From  the  electron  densities  alone,  measured 
by  the  spectrometer  in  the  discharge  and  extrapo¬ 
lated  from  the  radiometer  data,  there  appear  to  be 
more  electrons  [  resent  in  the  discharge  than  would 
be  predicted  by  the  Saha  equation;  e.g.,  an  argon 
plasma  under  total  thermodynamic  equilibrium 
(TE)  conditions  at  9000  K.  has  an  electron  density 
of  1.8  to  5.1  x  1015  e/cm3  for  pressures  of  100-700 
Torr.  These  values  are  considerably  less  than  the 
spectrometrically  measured  values  (Fig.  8).  Thus 
the  plasma  is  not  in  total  TE  over  the  range  of  dis¬ 
charge  pressures  considered.  However,  a  two- 
temperature  “bottleneck  level”  plasma  model  s?us- 
factorily  explains  the  experimental  observations. 
This  concept  will  be  elaborated  in  a  later  section. 

(b)  Trace  quantities  of  nitrogen  gas  (0.042  and 
0.084  mole  %)  were  injected  into  the  argon  to  ob¬ 
tain  N2  rotational  temperatures  in  the  channel  by 
using  the  rotational  lines  of  N2  and  N2+.  Specifi¬ 
cally,  the  rotational  lines  of  the  R  branches  of  the 
C3nu  -  B3rL  0-1  second  positive  band  of  N2 
were  observed  spectrometrically  in  the  channel 
when  the  power  was  54  kW  and  chamber  pressure 
139  Torr.  Uninverted  scans  of  the  unresolved 
triplets  in  the  R  branches  (J  =  32  to  J  =  45)  were 
treated  as  described  by  Robinson22.,  The  resulting 


rotational  temperature  was  calculated  to  be  1950K. 

In  addition,  the  rotational  lines  of  the  R  branch 
of  the  B^2U+  -  X3Sg+  0-0  first  negative  band  of 
N2+  were  observed  spe  trometrically  in  the  chamber 
at  48  kW  power  and  1 24  Torr  chamber  pressure. 
Nitrogen  gas  was  injected  into  the  argon  at  0.57 
mole  %.  An  uninverted  scan  was  made  and  treated 
as  described  in  Appendix  B.  Most  of  the  lines  from 
J  =  2  to  J  =  48  .vere  visible  but  weak  when  com¬ 
pared  with  the  intense  continuum  of  argon.  The 
rotational  temperatures  were  estimated  to  be  at 
least  2100  K. 

(c)  Calculations  of  the  gas  temperature  based 
on  gas  dynamic  considerations  indicate  that  the 
approximate  spatially-averaged  value  of  the  temper¬ 
ature  in  the  discharge  chamber  varies  between 
2900-4400  K  for  chamber  pressures  between  50-600 
Torr  (Fig.  13).  For  this  calculation  an  isentropic 
flow  from  the  discharge  chamber  was  assumed  with 
negligible  wall  loss.  Also,  the  flow  velocity  was 
assumed  to  be  sonic  at  the  discharge  chamber  exit 
with  the  mass  of  the  gas  being  carried  by  the  heavy 
gas  particles. 

(d)  Measurements  with  t  angsten-rhenium  ther¬ 
mocouples  in  the  channel  for  mass  flows  of  0. 1-0.36 
g/sec  and  42-ftO  kW  input  power  indicated  that  the 
gas  temperatures  ranged  between  1900-3320  K. 

(Fig.  1 4)^,;  The  mean  gas  temperature  calculated 
for  the  discbjfge  and  the  values  extrapolated  from 
the  thermocouple  data,  i.e.,  2500-4400  K,  are  in 
agree r  ient  if  isentropic  expansion  is  assumed.  The 
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Fig.  13  Estimate  of  peak  discharge  chamber  gas 

temperature  (calculated  from  gas  dynamics) 


Fig.  14  Gas  temperatures  in  reaction  channel 


thermocouple  measurements  as  well  as  the  calcu¬ 
lated  mean  temperatures  noted  in  c,  are  a  conserva¬ 
tive  lower  limit  to  the  true  centerline  temperatures. 
Comparison  of  all  of  these  values  to  the  spectrom- 
etrically  measured  excitation  and  rotational  tem¬ 
peratures  in  the  discharge  chamber  and  the  reaction 
channel  shows  that  TE  probably  does  not  exist  in 
either  place. 

(e)  As  was  noted  in  Section  3,  the  optical  and 
interferometer  data  from  the  channel  yield  closely 
correlated  values  for  the  electron  density  in  the 
channel  (Figs.  1  la  - 1  Id).  The  channel  electron 
density  varies  from  0.6  to  1.75  x  10*  ^  e/cm^  for 
the  range  of  discharge  conditions  considered. 

(f)  Optical  measurements  in  the  discharge 
chamber  at  54  kW  indicate  that  the  excitation  tem¬ 
perature  decreases  rapidly  from  10,000  K  between 
100-200  Torr  and  then  attains  a  steady  value 
around  8000  K  for  the  higher  pressures  (  >  300 
Torr).  Similar  behavior  has  been  noted  at  the 


other  power  levels.  The  measured  values  of  the 
excitation  temperature  are  in  agreement  with  those 
found  by  Scholz  and  Anderson2-*  at  75  and  380 
Torr  for  a  flowing  argon  rf  plasma  exhibiting  LTE 
(Fig.  15).  Their  discharge  was  operated  at  lower 
mass  flow  rates  and  power  levels  (see  Table  3).  The 
excitation  temperatures  here  were  obtained  from 
a  Boltzmann  graph  using  nine  line  intensities; 
Scholz  and  Anderson  used  four  lines.  Their  equilib¬ 
rium  temperature  was  obtained  from  the  absolute 
continuum  intensity  a‘  4259  A. 

(g)  The  excitation  temperatures  in  the  reaction 
channel  are  found  to  be  a  weakly  decreasing  func¬ 
tion  of  pressure  (for  alt  levels  of  rf  excitation)  with 
final  values  averaging  4500  K  (Table  2).  As  in  the 
discharge  chamber,  the  ratio  of  excitation  temper¬ 
ature  to  gas  temperature  (obtained  from  a  thermo¬ 
couple)  is  approximately  two  and  remains  essen¬ 
tially  invaiiant  for  the  entire  range  of  chamber 
pressures  and  rf  excitation  levels,  again  indicating 
that  the  plasma  is  in  a  two  temperature  state. 
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4.1  Two-temperature  bottleneck- level 
plasma  model 

The  foregoing  experimental  observations  and  analyt¬ 
ical  considerations  suggest  that  a  consistent  descrip¬ 
tion  of  the  plasma  both  in  the  discharge  and  the 
reaction  channel  requires  the  two-temperature 
bottleneck-level  plasma  model  suggested  in  a.  This 
model,  specifies  explicitly  that  excited  levels  above 
some  critical  energy  level  are  populated  in  equilib¬ 
rium  with  the  free  electrons  (at  temperature  Te). 

The  gas  is  at  lower  temperature,  Tg.  The  critical 
level  can  be  found  using  the  standard  criteria  for 
LTE  stated  by  Griem2*  for  hydrogen  and  hydrogen¬ 
like  atoms  or  ions.  Since  the  energy  level  diagram  of 
argon  is  similar  to  that  of  hydrogen,  these  criteria 
are  applicable  to  the  present  discussion. 

First,  the  LTE  criterion  for  the  upper  excited 
states  is  applied  (Griem,  Eq.  6-55).  Using  the  spec- 
trometric  results  for  the  discharge  chamber  of 
Ne  *  1  x  1016  e/cm3  and  T  2?  8000  K  the 
critical  level  for  LTE  is  shown  to  be  at  the  principal 
quantum  number  n  2  1.5.  Since  this  value  is  in¬ 
admissible  as  a  principal  quantum  number  but  the 
next  higher  integer  is,  it  is  concluded  that  all 
hydrogen-hke  levels  n  >  2  are  in  LTE,  Applying 


this  result  to  argon  requires  that  all  levels  from  the 
first  excited  state  on  are  in  LTE.  Similarly  applying 
the  Griem  criterion  to  the  channel  region  for 
Nc  s  3  x  1 0 1  ^  e/cm3  and  T  s  4500  K  gives  the 
critical  level  as  n  s  2.0  Thus,  the  first  excited  level 
of  argon,  'ing  somewhat  higher  than  that  for 
hydrogen, »  in  LTE  with  everything  above  it. 

For  the  ground  state  the  existence  of  LTE  must 
be  determined  separately  because  of  greater  devia¬ 
tions  from  hydrogenic  behavior.  For  these  cases  a 
modified  criterion  for  LTE  may  be  used  (Griem, 

Eq.  6-60)  which  requires  that  the  ground  level  be 
populated  ten  times  more  often  by  collisions  than 
by  radiation.  Applying  this  relationship  indicates 
that  the  electron  density  is  not  sufficient  to  guaran¬ 
tee  an  LTE  population  for  the  argon  ground  level. 
This  observation  supports  the  conclusions  obtained 
from  the  previous  LTE  criterion  above. The  fact 
that  the  Griem  criterion  is  satisfied  for  n  >  2  in 
both  the  discharge  chamber  and  the  reaction  channel 
“justifies”  the  interpretation  of  the  slope  of  the 
Boltzmann  graph  as  an  excitation  temperature. 
Finally,  the  hypothesis  of  LTE  is  also  supported  by 
the  close  agreement  of  the  electron  and  the  excita¬ 
tion  temperatures  as  obtained  from  the  spectrom- 
etric  and  radiometric  data  and  the  emissivity 
calculations. 
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5  Summary 


The  plasma  in  the  rf  discharge  and  the  sonic  after¬ 
glow  regions  appears  to  be  in  a  two-temperature 
state.  Comparison  of  the  electron  densities  obtained 
from  the  Kramers- Unsold  formulation  with  those 
obtained  f/om  an  equilibrium  calculation  indicates 
that  the  excitation  and  the  gas  (heavy  particle)  tem¬ 
peratures  are  not  equal.  However,  the  excitation  and 
electron  temperatures  are  approximately  equal  in 
the  discharge  itself,  averaging  8000  K  over  the  pres¬ 
sures  and  rf  power  levels  considered.  The  gas  tem¬ 
peratures  calculated  from  gas  dynamic  considerations 
and  extrapolated  from  thermocouple  data  average 
4400  K.  In  the  sonic  afterglow  region  the  corre¬ 
sponding  values  of  the  excitation  and  he?  vy  particle 
temperatures  are  almost  50%  less  than  those  meas¬ 
ured  in  the  discharge  chamber.  Independent  micro¬ 


wave  interferometer  measurements  of  the  electron 
density  confirm  the  Kramers- Unsold  values  obtained 
in  the  reaction  channel.  Optical  thinness  was  estab¬ 
lished  in  both  the  discharge  and  the  sonic  afterglow 

regions.  Using  the  standard  tests  for  LTE,  the  plasma 
is  shown  in  both  regions  to  be  in  LTE  in  all  but  the 
ground  level. 
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7  Appendix  A  -  Emissivity  calculation 


To  relate  the  effective  blackbody  radiation  temper¬ 
ature  of  i  plasma  source  to  the  electron  temperature, 
the  plasma  emissivity  must  be  known.  The  pertinen. 
expressions  for  the  plasma  emissivity  are  derived 
below. 

If  full  reciprocity  exists  in  the  propagation  char¬ 
acteristics  of  the  medium,  then  the  emissivity  and 
absorptivity  may  be  used  inte'  .  hangeably  as  a  con¬ 
sequence  of  KirchhofPs  law  of  radiation.  Mathe¬ 
matically  the  emissivity  (or  absorptivity)  may  be 
stated  in  the  usual  notation  as:7 


(  '  A  = 


/v  Re  ( T  •  E*)  dv 
/  Re  (H0  x  Hp  ds 


where  V  and  S  are  the  volume  and  projected  area  of 
the  sample  being  illuminated  by  an  electromagnetic 


test  wave,  composed  of  electric  and  magnetic  field 
intensities,  EQ  and  HQ,  respectively.  The  quantities 
J  and  E  are  the  current  density  and  electric  field 
intensity  within  the  sample.  Alternately  the  emis¬ 
sivity  may  be  formulated  as 

e  =  A  =  1-R— T  (A-2) 


where  T  and  R  are  the  power  transmission  and 
reflection  coefficients,  respectively.  Specific  ex¬ 
pressions  for  R  and  T  will  be  derived  for  a  slab 
plasma. 

For  the  case  where  the  plasma  can  be  considered 
a  plane,  effectively  homogeneous  slab,  a  plane  wave 
analysis7 ,  assuming  incoherence  in  the  reflected 
waves  within  the  slab,  yields  a  reflection  coefficient, 


\  /  /Kj  ' 

-^r3j  cosk2«+j  {/ 


+  \nr3)  cos  k2^  +  j 


sin  k-,5 


K3  i?o  V 

,_^+VT2)  stak2e 


\T&2  =  ^2  “7X2;  V"K"3  =  ^3  ~->*3;  k2  ~  TT  V^2  . 
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is  the  wave  impedance  of  the  medium  outside  of 
layers  2  and  3,  tjQ  =  377  SI  (the  free  space  imped¬ 
ance),  and  K2  and  K3  are  the  complex  dielectric  con¬ 
stants  of  layers  2  and  3,  representing  in  the  present 
analysis  the  water  layer  in  the  discharge  jacket  and 
the  plasma,  respectively.  The  power  reflection 
coefficient  is  defined  by  R  =  |pj^  where  p  is  given 
in  Eq.  (A-3). 


I  [-0.35  +  0.775  (10— 4)  Te]  10"16; 
for  Te  >  104K 
[0.39  -  0.551  (10— 4)  Te 
+  0.595  (10“8)  Te2  ]  10  16; 
for  Te  <  104K.  (A-6) 


Measurements  indicate  that  the  plasma  attenua¬ 
tion  in  the  present  case  exceeds  50  dB;  thus  TsO. 
The  defining  equations  for  the  plasma  (layer  3)  are: 


**3 

*3 


(A-4) 


+ 


1/2 


An  additional  relationship  can  be  derived  for  Ne  in 
terms  of  the  discharge  chamber  pressure  and  mass 
flow  rate  giving 


0.962  x  1019P- 


Ne  = 


0.7545  =  10_4 /£tV 


M, 


M0 

- -  +  Te 

1  +  p  e 


(A-7) 


Five  assumptions  were  made  in  deriving  Eq.(A-7): 

1 )  isentropic  flow  from  the  discharge  chamber,  2) 
sonic  flow  at  the  exit,  3)  applicability  of  the  perfect 
gas  law,  4)  mass  of  the  gas  carried  by  the  heavy 
particles,  an  i  5)  no  significant  temperature  gradient 
in  the  chamber. 


The  collision  frequency  for  electron-ion  inter¬ 
actions  (Coulomb  interactions),  vQV  is  given  using 
the  Appleton-Bray  formulation2  as 


ei 


°  <2kTe)3/2 


fin 


k3Te3  1 
e 

(A-5) 


The  foregoing  set  of  equations  defining  the  plasma 
emissivity  for  the  slab  case  may  be  rearranged  for 
computer  calculations  so  that  the  electron  temper¬ 
ature  and  the  percentage  of  ionization  may  be 
retained  as  parameters  in  the  calculations.  A  typical 
result  from  these  calculations  is  given  in  Fig.  1 2, 
indicating  reaconable  agreement  with  the  experi¬ 
mental  data  from  the  radiometer  and  spectrometer. 


(as^u.-iing  that  Nj  =  Ne).  The  collision  frequency 
for  the  less  dominant  electron-neutral  interactions 
is  given  by 


v 


ea 


8kTe\  */ 2 

it  mey 


where  the  collision  cross  section,  oea,  is  given  as3 
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Symbols  used  in  Appendix  A: 


N| 

Ne 

o)p 

f 

V 

Te 

Py 

M 

P 

M0 

C 

c,e,me,k 


=  ion  density  (ions/cm-*) 

=  electron  density  (e/cm^) 

=  plasma  frequency  (<op  =  5.63  x  1 04\/N^) 
=  frequency  of  the  EM  wave  (Hz) 

=  »ea  +  i>ej  =  effective  collision  frequency 
=  electron  temperature  (K) 

=  discharge  chamber  pressure  (Torr) 

=  mass  flow  rate  (g/sec) 

=  fraction  of  ionization 
=  atomic  weight  (g) 

=  thickness  of  layer  2 
=  velocity  of  light  (cm/sec);  electron 
charge  (e.s.u.),  electron  mass  (g), 
Boltzmann  constant 
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8  Appendix  B  -  Determination  of  rotational  temperatures* 


The  electron  translational  temperature,  in  laboratory 
plasmas,  usually  determines  the  excitation  tempera¬ 
ture  of  the  heavy  particle  energy  levels.  It  is  of 
interest  to  learn  that  the  heavy  particle  translational 
temperature  can  determine  the  rotational  tempera¬ 
ture  of  a  molecular  species.  Thence  molecular  spec¬ 
troscopy  might  be  a  tool  for  finding  this  otherwise 
evasive  property.  For  this  reason,  the  present  deter¬ 
mination  of  a  molecular  rotational  temperature 
was  undertaken. 

The  plasma  jet  which  produced  these  spectra 
was  being  operated  in  the  laminar,  subsonic  mode, 
and  there  was  reason  to  hope  that  local  thermo¬ 
dynamic  equilibrium  existed.  The  0-0,  1-1 ,  and  2-2 
bands  of  N2  (region  3700-3900  A)  were  well  devel¬ 
oped  and  at  least  the  0-1  band  (region  4200-4300  A) 
was  visible.  All  of  these  bands  belong  to  the 
system,  which  is  called  a  “negative  system”  (in  con¬ 
trast  to  the  “positive  systems”  which  are  emitted 
by  neutral  N2)  because  of  where  it  appears  in  a 
common  discharge  tube.  The  appearance  of  molec¬ 
ular  nitrogen  is  due  to  the  seeding  of  the  plasma 
jet  operating  gas,  argon,  by  a  50%  (by  weight) 
aqueous  solution  of  Mn  (N03>2  in  such  an  amount 
that  there  were  1 2.9  nitrogen  atoms  per  million 
argon  atoms.  The  spectrometer  saw  light  from  a 
diameter  through  the  plume  1-in.  downstream  of 
the  exhaust  plume. 


8.1  Theory 


The  relative  intensities  of  entire  bands  emitted  by 
molecules  in  an  optically  thin  gas  in  local  thermo¬ 
dynamic  equilibrium  (LTE)  are  proportional  to  the 
populations  of  the  upper  rotational  levels,  given 
by  the  equation: 

n  (v')  =  n(0)e— )hc/kT,  (B-l) 


(Ref.  1 )  in  which 

G(v')  =  coe'  (v'  +  1/2)  +  we'  xe'(v'  +  1/2)2 
+  C0e'  ye'  (v'  +  1/2)3 

and 


(B-2) 


v'  =  vibrational  quantum  number  of  upper  levels, 
T  =  vibrational  excitation  temperature. 


coA',  c oj  x  '  00'  y  '  are  numbers  tabulated  by 
Herzberg  (Ref..  1). 


The  relative  intensities  of  the  rotational  lines  in 
an  optically  thin  gas  in  LTE  within  a  given  band 
and  branch  are  given  by 


I  =  C  —  (j'  +  r+De-Bv'J'H'+Ohc/kT 


<B-3) 


*Material  in  this  appendix  is  taken  from  an  internal 
McDonnell  Douglas  Report  entitled  "A  Method  of 
Obtaining  Rotational  Excitational  Temperature 
from  Spectroscopic  Observation  of  a  Plasma"  by 
R.J.  Checkley  and R.A.  Hefferlin (Dec.  1967). 


The  measurements  in  this  appendix  were  made  in  an  arc 
plasma  source,  that  was  entirely  different  from  the  rf 
discharge  described  earlier.  The  method  of  analysis  is  the 
same  as  described  in  Section  4. 
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for  which 

T  =  rotational  excitation  temperature, 

V  -  rotational  quantum  number  of  the  upper 

level, 

I"  =  rotational  quantum  number  of  the  lower 
level. 


The  selection  rules  are  such  that  for  a  P  branch 
(J'  +  J"  +  1)  =  (J"  -  1)  +  (J"  +  1)  =  21",  (B-4) 

and  for  an  R  branch 


(J'  +  J"  +  1)  =  (J"  +  1)  +  (J"  +  1)  =  2J"+2 

(B-5) 

Also 

By'  =  Be  —  are  (v'  +  1/2),  (B-6) 

where  v'  is  the  vibrational  quantum  number  of  all 
the  upper  rotational  levels  and  Be  and  Og  are  molecu¬ 
lar  constants  again  tabulated  by  Herzberg. 

Algebraic  manipulation  of  Eq.  (B-3)  yields  (base  10) 


IX4  By'  J'  (J'  +  1)  he 

108  (F+rTT)  _,°8C~  23  kT 

(B-7) 

If  the  entire  left  hand  side  is  plotted  on  a  y  axis,  and 
J'  (J'  +  1)  is  plotted  on  an  x  axis,  the  points  plotted 
should  fall  along  a  straight  line  with 

Bv'  he 

Slope  =  — -  (B-8) 

p  2.3  kT 

Solution  of  this  equation  then  gives  the  rotational 
excitation  temperature  (correction  must  be  made 
for  the  alternation  of  intensities  which  results  from 
the  influence  of  nuclear  spin;  otherwise  two  parallel 
lines  will  appear  on  the  temperature  graph  due  to 
the  fact  that  every  other  line  is  twice  as  bright  as 
its  neighbors.) 


In  the  particular  case  at  hand,  the  quantum 
number  J  "  is  replaced,  by  molecular  spectroscopists, 
with  the  quantum  number  K  ",  or,  identically,  with 
K.  The  quantum  number  J '  could  likewise  be  re¬ 
placed  by  K;  however  it  is  more  common  to  use 
the  selection  rules,  Eqs.  (B-4)  and  (B-5),  to  express 
J  'in  terms  of  K2,s.  We  retain  the  J  'and  J  "notation, 
however,  in  this  report. 

The  procedure  of  obtaining  the  rotational  tem¬ 
perature  presupposes  that  one  can  identify  the 
rotational  quantum  numbers  of  the  individual  lines. 
The  difficult  job  of  accomplishing  this  identification 
will  be  described  in  the  following  section,  as  will 
the  explanation  of  a  correction  applied  due  to 
blending  of  lines. 


8.2  Procedure 


At  low  temperatures  (e.g.,  room  temperature)  the 
(absorption)  spectrum  of  the  negative  B2  X-X2  £ 
bands  of  the  molecular  ion  exhibits  clear  sepa¬ 
ration  of  the  scpa.ute  bands,  thus  leaving  only  the 
problem  of  separating  lines  from  the  various 
branches  of  one  band  from  one  another/ 

At  higher  temperatures  the  emission  lines  of  the 
different  bands  lie  among  each  other,  causing  a 
much  more  difficult  problem.4  In  the  spectrum 
under  study,  the  lines  from  various  bands  were 
overlapped  to  the  extent  that  attempts  to  assign 
rotational  quantum  numbers,  branch  numbers,  and 
vibration  ll  quantum  numbers  by  comparison  with 
a  published  spectrum4  were  completely  unsuccess¬ 
ful. 

A  model  of  the  spectrum  was  therefore  con¬ 
structed  from  basic  molecular  data  and  from 
assumed  source  conditions.  Wavelengths  of  the 
first  1 50  lines  from  the  P  and  R  branches  of  the 
0-0,  1-1,  2-2,  and  3-3  bands  were  calculated  from 
molecular  constants  given  by  Herzberg.  Relative 
intensities  of  lines  in  selected  wavelength  regions 
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were  calculated  by  multiplying  (a)  the  exponential 
term*  of  Eq.  (B-l ),  (b)  the  parenthesis  term  and 
the  exponential  term*  of  Eq.  (B-3),  and  (c)  the 
intensity  alternating  factor  due  to  nuclear  spin.  A 
histogram  was  constructed  by  summing  all  these 
intensities  in  0.8  A  increments  of  wavelength.'*' 

This  histogram  strikingly  resembled  the  strip  chart 
spectra.  Those  few  of  the  histogram  peaks  which 
were  due  primarily  to  a  angle  line*of  a  single  branch 
of  one  of  the  bands  then  allowed  the  assignment  of 
rotational  quantum  numbers  to  the  corresponding 
spectrum  features.  AH  of  these  features  were  in  the 
0-0  band.  The  assignment  of  rotational  quantum 
numbers  was  checked  by  plotting  a  Fortrat  parabola 
(quantum  number  vs  wave  number)7 ;  the  wave 
number  minimum  of  the  parabola  was  found 
algebraically  to  fall  at  the  0-0  band  head.  These 
considerations  also  allowed  assignment  of  rotational 
quantum  numbers  in  the  R  branch  of  the  0-0  band 
to  a  very  regular  series  of  lines  far  to  the  blue  side 
of  all  the  band  heads  and  blended  lines  (J  "  =  70  and 
beyond).  It  was  next  necessary  to  measure  the  inten¬ 
sities  of  as  many  lines  as  possible.  There  would  not 
be  many  because  the  number  of  spectrum  features 
due  primarily  to  a  single  line  -  according  to  com¬ 
parison  with  the  histogram  -  was  small.  And  bending 
was  to  be  expected  even  in  the  features  now  identi¬ 
fied  with  J"  >  70.4 

The  intensities,  I,  were  obtained  by  multiplying 
the  heights  of  the  isolated  features  (top  of  strip 
chart  =  100)  by  their  widths  (full  width  at  half 
maximum  in  mm)  and  by  the  spectrograph  calibra¬ 
tion  obtained  with  a  tungsten  standard  lamp.  These 
intensities  were  then  dealt  with  in  two  ways: 

1.  They  were  multiplied  by  X4,  and  by  2  if  J  " 
was  odd;  they  were  divided  by  2  (J  "  +  I )  if  in  the 
R  branch  and  by  2J  "if  in  the  P  branch; and  the  com¬ 


*It  was  assumed  that  T  was  12,500  K 
+Our  instrument  bandpass 

**Due  to  the  low  temperature  obtained  in  the  first  method 
the  correction  factors  of  Ref.  (4)  were  linearly  extrapolated 
to  4000  K.  Linear  interpolation  was  considered  satisfactory 
as  far  as  the  quantum  numbers,  J  ”,  were  concerned. 


mon  logarithm  was  taken.  The  results  of  these  opera¬ 
tions  were  plotted  against  J’(J+ 1),  as  required  by  Eq. 
(B-6).  A  typical  calculation  appears  in  Table  B-l. 

The  plot  is  shown  in  Fig  B-l ;  the  slope  of  the 
selected  straight  line  is  -3.36  x  10  . 

2.  They  were  multiplied  by  X*  and  by  the  correc¬ 
tion  factors**  given  in  (Ref.  4);  they  were  divided  by 
2  (J  "  +  1 )  if  in  the  R  branch  and  by  2  J  ”  if  in  the 
P  branch;  the  common  logarithm  was  taken.  The 
results  of  these  operations  were  plotted  against 
J  (J  +  1 ),  as  required  by  Eq.  (B-6).  A  typical  calcula¬ 
tion  appears  in  Table  B-l. The  plot  is  shown  in 
Fig.  B-2;  the  slope  of  the  selected  straight  line  is 
-3.84  x  10~4. 


8.3  Results 

The  use  of  Eq.  (B-7)  made  it  possible  to  find  the 
rotational  excitational  temperatures  for  the  two 
methods  once  the  slopes  of  the  Boltzmann  plots 
(Figs.  B-l  and  B-2)  were  known.  Respectively, 
these  temperatures  were  4300  K  and  3750  K. 

Due  to  the  subjective  nature  of  the  selection  of 
the  lines  in  Figs.  B-l  and  B-2,  these  numbers  are 
considered  to  mean  that  the  correct  value  is  4000  K 
+  300  K.  This  rotational  excitational  temperature 
is,  of  course,  an  average  through  the  plume  of  the 
plasma  jet. 


8.4  Conclusions 

The  rather  large  scatter  of  the  points  in  Figs.  B-l 
and  B-2  is  probably  due  (a)  to  fluctuations  in  the 
source  and  (b)  to  the  presence  of  one  or  two 
atomic  spectrum  lines. 

The  histogram,  which  makes  quantum  number 
identification  possible,  does  not  depend  upon 
having  presupposed  the  correct  temperature  in  the 
exponential  terms  of  Eq.;  (B-l )  and  (B-3). 
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Table  B-l 

Typical  calculations 

(Using  0-0  bard,  rotational  line  J"=  79,  in  R  branch) 


Numerical 

calculation 

Quantity  or  operation 

Using  Eq.  (B-& 

Using  eon  action 
factor  of  Ref.  4 

(1)  Line  recording  area  (top  of  chart  =  100; 
FWHM  in  mm) 

13.5 

13.5 

(2)  Instrument  response  function  (relative) 

2.25 

2.25 

(3)  Wavelength  (in  A) 

3,725 

3,725 

(4)  X4 

1.92  x  1014 

1.92  x  1014 

(5)  Multiple  by  2  if  J  "  is  odd 

(6)  Correction  factor  of  (Ref.  4) 

2 

0.88 

(7)  Divide  hy  2(J "  +  1 )  if  in  R  branch,  by 

2J  "  if  in  P  branch 

160 

160 

(8)  Product  of  (1),  (2),  (4).  (5).  or  (6), 
divided  by  (7) 

7.2  x  1013 

3.2  x  1012 

(9)  Logarithm  (base  10),  which  is  plotted  on 

13.862 

12.505 

"y"  axis  of  Fig.  B-1  or  Fig.  B-2 

do)  j'  • 

80 

80 

(11)  J(J'+  1).  which  is  plotted  on  "x"  axis  of 
Fig.  B-1  or  Fig.  B-2 

6,480 

6,480 

Although  two  methods  were  used  to  reduce  the 
raw  data,  one  of  which  assumed  that  even  those 
features  finally  isolated  were  still  blended,  these 
two  methods  gave  substantially  the  same  final 
result. 
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